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ABSTRACT 



Context. This paper is a follow-up on the vast effort to collect radial velocity data for stars belonging to the Hipparcos survey. 
Aims. We aim at extending the orbital data available for binaries with M giant primaries. The data presented in this paper will be used 
in the companion papers of this series to (i) derive the binary frequency among M giants and compare it to that of K giants (Paper II) 
and (ii) analyse the eccentricity - period diagram and the mass-function distribution (Paper III). 

Methods. Keplerian solutions are fitted to radial-velocity data. However, for several stars, no satisfactory solution could be found, 
even though the radial- velocity standard deviation is greater than the instrumental error, because M giants suffer from intrinsic radial- 
velocity variations due to pulsations. We show that these intrinsic radial-velocity variations can be linked with both the average 
spectral-line width and the photometric variability. 

Results. We present an extensive collection of spectroscopic orbits for M giants with 12 new orbits, plus 17 from the literature. 
On top of these, 1 preliminary orbit yielded an approximate value for the eccentricity and the orbital period. Moreover, to illustrate 
how the large radial- velocity jitter present in Mira and semi-regular variables may easily be confused with orbital variations, we also 
present examples of pseudo-orbital variations (in S UMa, X Cnc, and possibly in HD 1 15521, a former IAU radial-velocity standard). 
Because of this difficulty, M giants involving Mira variables were excluded from our monitored sample. We finally show that the 
majority of M giants detected as X-ray sources are actually binaries. 

Conclusions. The data presented in this paper considerably increase the orbital data set for M giants, and will allow us to conduct a 
detailed analysis of the eccentricity - period diagram in a companion paper (Paper III). 

Key words, binaries: spectroscopic - stars: late-type 



1. Introduction 

When a mass-losing M giant is present in a binary system, 
the interaction of the wind of the giant with the compan- 
ion gives rise to photometric activity or spectroscopic symbio- 
sis (like in symbiotic stars, in VV Cephei-like systems, or in 
low-mass X-ray binaries like V2116 Oph) that make the sys- 
tem very conspicuous even far away. However, if the M giant 
does not lose large amounts of mass, the binary nature of the 
star will not be as conspicuous. Combined with M giants of- 
ten exhibiting pulsations that cause intrinsic velocity variations, 
thereby confusing the search for orbital variations, it explains 
why so few spectroscopic binaries involving M giants are known 
so far. Inde ed, the Ninth Catalo gue of Spectroscopic Binary 
Orbits (Sr9 : iPourbaix et al] |2004) contains 2746 entries (query 
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in March 2007), among which only 32 systems involve M giants, 
and yet 21 of these are either well-known symbiotic systems or 
VV-Cephei-like systems. 

Carquillat and collaborators have devoted a series of pa- 
pers to spectroscopic binaries of spectral types F to M, but 
only 2 binaries with primaries of spectral type Mi ll were stud- 
ied (ICarquillat & Ginestetlll996l:lPrieuret al.ll2006l) . Previously, 
IStephensonl (1 19671) provided a list of 7 systems wi th composite 
spectr a involving an M star. But it is the paper by Hinkl e et~ai1 
(2002) that, to the best of our knowledge, has so far provided 
the most extensive list of spectroscopic binaries involving non- 
symbiotic M giants. 

With the present paper, we start a series devoted to a de- 
tailed study of the properties of spectroscopic binaries involv- 
ing an Mill primary. The number of such binaries with known 
orbital elements has nearly doubled, as a result of our observ- 
ing campaign of an extensive sample of M giants, drawn from 
the Hipparcos Catalogue, for which CORA VEL radial veloci - 
ties have been obtained in a systematic way dUdrv et al.l fl997). 
The main driver behind this large database lies, of course, with 
the stellar kinematics in our Galaxy. And indeed the kinemati- 
cal properti es of the presen t sample of M giants have been fully 
analysed bv lFamaev et al.1 d2005l) and lFamaev et all (120081) . But 
this large data set may also be used to search for binaries. 

The present paper presents the radial-velocity data (Sect. 2) 
and orbital elements (Sect. 4) of the newly-discovered spectra- 
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scopic binaries, for which a satisfactory orbit could be obtained. 
It also discusses the intrinsic variations sometimes mimicking 
orbital variations (Sect. 3). The list of new orbital elements is 
complemented with an exhaustive list of orbital elements for 
non-symb iotic M giants drawn f rom the literature (Table [§). In 
Paper II (Frankows ki et al.ll2009l) . we will use the observational 
information gathered in this paper to derive the frequency of 
spectroscopic binaries among M giants, and compare it with that 
of K giants. Paper III (Jorissen et al. 2009) will then present an 
in-depth analysis of the eccentricity-period diagram for M gi- 
ants. 



2. Radial-velocity data 

The basic sample of M giants is drawn from the Hipparcos sur- 
vey stars (identified b y flag 'S' in field H68 of the Hipparcos 
catalogue; ESAI I19971) . These M giants were extracted from the 
Hipparcos survey stars on the basis of the spectral type pro- 
vided in the Hipparcos catalogue and of the absolute magni- 
tude Mjj p < 4 computed from the Hipparcos parallax and Hp 
magnitude. Mira stars or supergiants of luminosity class I (when 
explicitly mentioned in the spectral classification) were not in- 
cluded in the sample (notably because of the confusion that their 
envelope pulsation may cause on the radial-velocity variations). 
This first sample (defined here as sample I) contains 771 M gi- 
ants with declinations greater than —15° and corresponds to the 
sample of northern M giants analysed in Famaey et al. (2005), 
to which 65 M giants with declinations between and —15° 
have been added. Two radial-velocity measurements, spanning 
at least one yea r, have been obtain ed with the CORAVEL spec- 
trovelocimeter (Ba ranne et al.ll 1 979b for all stars of sample I, as 
part of a monitoring pro gramme targeting all Hipparcos survey 
stars later than about F (|Udrv et aklll997l) . Note that 22 objects 
of the sample of lFamaev et alj d2005l) are not present in sample I 
because they had only one radial velocity measurement, making 
them unsuitable forbinarity analysis. Famaey's sample was also 
screened for other irregularities, such as wrong Hipparcos spec- 
tral type or mistaken identity^. Extrinsic S stars were removed 
from the sample as well. 

Subsamples of sample I have been subject to more intensive 
observing campaigns. Every third star from this sample has re- 
ceived a denser coverage, with (mostly) 4 instead of 2 measure- 
ments (7 stars received only 3 measurements), to achieve a bet- 
ter binary detection rate: this subsample of sample I is defined as 
sample II. Figure Q] displays the histograms of N mcSiS (the num- 
ber of measurements per star), < e > (average uncertainty of one 
measurement), and At (time spanned by the measurements of a 
given star) for the 254 stars of sample II. 

Furthermore, 35 stars from sample II, suspected of being bi- 
naries (i.e., with a radial-velocity standard deviation a(V r ) > 
1 km s _1 at the end of the observing campaign of s ample II) have 
been monitored with the ELODIE spectrograph (Ba ranne et alj 
1996) at the Haute-Provence Observatory to derive their orbital 
elements. They make up sample III. At the end of the ELODIE 
monitoring, a late measurement was obtained for 157 stars from 
sample II with cr(V r ) < 1 km s _1 (this selection being mainly 
based on right ascension), in order to detect binaries with very 



1 The only mistake found in Famaey et al] d2005t ) was the star 
HIP 26247 = RR Cam = BD +72° 275, which was wrongly as- 
signed radial velocity measurements from the binary star J275 in the 
Hyades cluster. The former star should therefore be discarded from the 
CORAVEL sample. 
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Fig. 1. Histograms of the number of measurements per star 
(iV meas ), the average uncertainty of one measurement (< e >), 
and the time spanned by the measurements of a given star (Ai), 
for sample II (see Tables Q] and |2). 

Table 2. The first five lines of the list of stars from sample II. 
The full table is only available electronically at the CDS. 



HD JV tot iV EL o V T ff(VV) e(VV) W fte^ 
(kms -1 ) (kms -1 ) (kms^ 1 ) (km s" 1 ) 



1013 


5 


1 


-46.67 


0.18 


0.26 


2.80 


1255 


5 


1 


14.57 


0.32 


0.27 


3.95 


2313 


5 


1 


24.15 


0.18 


0.28 


3.15 


2411 


6 


2 


4.67 


1.06 


0.25 


4.58 


2637 


5 


1 


11.25 


0.19 


0.25 


1.33 



ORB 



a The average sigma of a Gaussian fitted to the cross-correlation profile, 
corrected for the instrumental width (7 km s _1 for CORAVEL at the 
Haute-Provence Observatory); 
b The flag in this column is the same as in Table 4. 



long orbital periods. These constitute sample IV. A summary of 
the properties of these four samples is presented in Table Q] 

The CORAVEL data were put o n the ELODIE r adial- 
velocity system to ensure homogeneity (lUdrv et al.l fl999b). The 
uncertainty < e > of one radial-velocity measurement is approx- 
imately 0.3 km s _1 for CORAVEL measurements, but is bet- 
ter for ELODIE measurements. With the ELODIE single-fibre 
mode used during th e present observing campaigns, it may be 
as good as 50 m s _1 dBaranne et al . 1996). It was not measured 
in a systematic way; however, to fix the ideas, an accuracy of 
0.2 km s _1 has been associated with the ELODIE measurements 
in the data files. 

The data for sample I (average velocity, radial-velocity stan- 
dard deviation, and binarity flag) may be found in Table A.lQ 

2 As a result of the confusion between HIP 26247 and a binary star 
(see footnote 1), the star HIP 26247 is erroneously flagged as a binary 
in that table. 
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Table 1. The four observing campaigns defining the four (sub)samples of Hipparcos M giants. 
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Sample 



III 



IV 



Selection 2 measurements in Famaey et al. (2005) 



iVstar 

iVmcas per star 
Telescope 

Spectrograph 
Time span 



771 

2 

Swiss 1-m (OHP) 
CORAVEL 
1991-1999 



4/3 of I allcr(V;) > Hems" 1 from II a(V r ) < 1 km s" 1 from II 



254 
> 4 a 
Swiss 1-m (OHP) 
CORAVEL 
1993-1995 







35 
> 5 
1.93-m (OHP) 
ELODIE 
2000-2004 



* < a < 16.5" 
138 
7V(II)+1 
1.93-m (OHP) 
ELODIE 
2001-2004 



Seven stars have N n 



= 3. 



Table 3. The individual radial-velocity measurements for the 35 
stars of sample III. The full table is available electronically from 
the CDS, Strasbourg. 



HD 


JD-2400000 


Vr 


e(Vr) 


Inst" 


2411 


48876.599 


4.97 


0.32 


COR 


2411 


49223.556 


6.05 


0.33 


COR 


2411 


49668.360 


6.42 


0.31 


COR 


2411 


50717.506 


4.74 


0.32 


COR 


2411 


52304.297 


3.56 


0.20 


ELO 


2411 


53047.246 


4.39 


0.20 


ELO 


4301 


48441.832 


6.36 


0.28 


COR 


4301 


49309.361 


7.08 


0.29 


COR 


4301 


49581.863 


6.96 


0.27 


COR 


4301 


49602.549 


6.58 


0.29 


COR 


4301 


53048.254 


4.50 


0.20 


ELO 



a A flag identifying the spectrograph: COR = CORAVEL; ELO = 
ELODIE; CAM = Cambridge spectrovelocimeter 

(columns 24-26) of lFamaev et alj d2005l) . A similar table, merg- 
ing CORAVEL and ELODIE data, is given in Table |2] for sam- 
ple II, which also provides the binarity diagnostics (according 
to the rules that will be specified in Sect. 13.11 ). Individual radial- 
velocity measurements for the 35 stars of sample III are given 
in Table [3] Roger Griffin has kindly provided us with supple- 
mentary measurements for HD 182190 and HD 220088 that al- 
lowed us to compute an orbit for these two stars. These measure- 
ments were performed with his Cambridge spectrovelocimeter 
and are listed in Table[3] with the label 'CAM'. To put them on 
the ELODIE system, an offset of —0.8 km s _1 has been applied. 

3. Binaries, intrinsic radial-velocity jitter, and 
pseudo-orbits caused by pulsation 

3.1. Binarity diagnostic 

The search for spectroscopic binaries (SBs) among M giants is 
made difficult by the bulk mass motion existing in the atmo- 
sphe res of these stars (all M giants being variable to some extent; 
e.g. [ Ever & Grenonlll997t iJorissen et al.lll997t ; rSoszvhsk i et al.l 
2004), since such motion triggers some intrinsic radial-velocity 
jitteiS (e.g., lUdrv et al.|[T99l iHinkle et al.ll2002l) . This intrin- 
sic jitter amounts to nearly 1.5 km s _1 in the coolest (non- 
Mira) M giants^ (Fig. |2J, so the detection of SBs cannot rely 

A term introduced bv lGunn&Griffid ( ll979l) in this context. 
4 Mavor et al] i 19841) have shown that stars located at the tip of the gi- 
ant branch in the globular cluster 47 Tuc have a jitter of ~ 1.25 km s - , 
while the jitter reduces to 0.27 km s _1 for 1 mag-fainter stars. This trend 
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Fig. 2. The radial-velocity standard deviation cto (Vr) as a func- 
tion of the average intrinsic line width Sb, for samples II, III, and 
IV (see Table [TJ- For non-binary stars (open squares), the max- 
imum radial-velocity jitter increases with Sb. The dashed line 
represents a rough division between binary and non-binary stars. 
How this line has been defined will be discussed in Paper II. 
Symbols are as indicated in the figure. 



solely on a \ 2 test comparing the radial-velocity standard de- 
viation cr(Vr) to the average instrumental error < e >. In the 
following, we denote by <To(V r ) the radial-velocity standard de- 
viation corrected for the average instrumental error < e >, i.e. 
ao(K) = (^ 2 (K)-<e> 2 ) 1/2 - 

In a first step, to flag a star as a binary in samples II, III, 
and IV, we therefore rely solely on the visual examination of the 
radial-velocity variations, and whether or not it is possible to ob- 
tain a meaningful orbital solution. The binarity diagnostic is then 
complemented by two criteria, described in detail in Sects. [37X7X1 
and |3.1.2l (i) the location of the star in a diagram <Jo(V r ) — Sb 
(Fig. |2), where Sb denotes the intrinsic width of spectral lines 
(Sect. [37T7iT ), and (ii) its location in a diagram tJo(V r ) — cr(Hp) 
(Fig. [3}, where <r(Hp) is the standard deviation of the Hipparcos 
magnitude (Sect. |3.l72l . The final binarity diagnostic for all stars 



goes further do wn the red-giant branch , as shown with much hig her ac- 
curacy levels bv lSetiawan et al.1 d2004h and lda Silva et al.1 d2006h . 
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Fig. 3. The radial-velocity standard deviation (corrected for the 
instrumental error) as a function of the standard deviation on the 
Hipparcos Hp magnitude. Symbols are the same as in Fig. [2] 
except that large circles identify stars detec ted as X-ray source s 
by the ROSAT All-Sky Survey, following iHiinsch et al.l (1 19981) . 
The dotted line represents the trend between radial-ve locity and 
photometric variability reported by iHinkle et al.l d 1 997b . This is 
not a dividing line like the one plotted in Fig. 2. 



from sample III (and for stars from sample IV, which are SB or 
suspected SB) is listed in Table|4]according to the following cat- 
egories: 

a "ORB" : a satisfactory orbit could be computed (see Fig. 4); 

b "ORB:" : the orbit is preliminary, because the number of data 
points or the time coverage are not large enough, or two dif- 
ferent solutions are possible (see Fig. 4); 

c "SB" : the star is a spectroscopic binary, but there are not 
enough data points to even compute a preliminary orbit (see 
Fig. 5); 

d "SB?" : the star is suspected to be a binary (see Fig. 6), even 
though it falls close to the dividing line in the ao(Vr)—Sb 
diagram (see Sect. |3.1.T| and Fig.|2|. 

e "SB / jitter?" : it is not clear whether the radial velocity vari- 
ations are of intrinsic or extrinsic nature (see Fig. 7): the star 
falls close to the dividing line in the ao(Vr)-Sb diagram 
(Fig.©. 

The radial-velocity variations (of intrinsic nature) of non-SB 
stars of sample III are plotted in Fig. 8. 

3.1 .1 . Intrinsic width of spectral lines 

We note that the parameter Sb, measuring the intrinsic width of 
spectral lines, provides useful guidance in distinguishing radial- 
velocity jitter from orbital motion. This parameter is defined 
as Sb = (s 2 — so 2 ) 1 / 2 , where s is the standard deviation of 
a Gaussian fitted to the CORAVEL cross-correlation dip (see 
Baranne et al. 1979) and sq the instrumental width (7 km s _1 
for CORAVEL at the Haute-Provence Observatory; for details 



Fig. 4. Radial-velocity data points for the spectroscopic bina- 
ries, along with their orbital solution. The solution marked as 
'ORB:' is not very well constrained and hence preliminary. All 
panels span a Vr range of 30 km s _1 . Error bars on the individ- 
ual measurements (of the order of 0.2 km s -1 ) are too small to 
be seen. 



see Uorissen & Mavorl 119881: IVan Eck & Jorissenl 120001 also 
Paper II). Indeed, Fig. [2] reveals that stars for which no orbital 
solution could be found (open squares) have radial-velocity stan- 
dard deviations that increase as Sb increases. In Paper II, we will 
show that Sb correlates closely with the stellar radius. We thus 
foresee that radial-velocity variations in stars falling well below 
the dashed line in Fig. [2] are most likely of intrinsic (rather than 
orbital) origin. Of course, binaries with an orbit seen almost face 
on, or with a long period, or insufficiently sampled, may also 
fall below the dashed line. Let us stress that the slope of this di- 
viding line in Fig. [2] can be influenced by both the number of 
radial-velocity data points for a given object and the total num- 
ber of objects in the sample (see Sect. 2.1 of Paper II). 

Examples of radial-velocity variations for such stars from 
sample EH, falling below the dividing line in Fig. [2j are plotted 
in Fig. [HJ and it is clear that no orbital solution can be found 
for them. On the other hand, stars along the dashed line, whose 
velocity variations are plotted in Fig. [7] could either be SB or 
exhibit radial-velocity jitter. These 10 stars are listed in Table|4]:. 

On the other hand, among stars from samples III and IV, 
for which 5 or more measurements have been made, 22 certain 
spectroscopic binaries have been found, as listed in Table|4^. Of 
those, only 3 (HD 89758, HD 108907, and HD 132813) were 
previously known to be spectroscopic binaries, and a new satis- 
factory orbit could be computed for 12 of them (with the help 
of supplementary data points from the Cambridge spectrove- 
locimeter for HD 182190 and HD 220088; see Tabled The 
orbital elements are presented in Sect. [4] The radial-velocity data 
points of the binary stars are displayed in Fig. [4] for those bina- 
ries with an orbit available (labelled in Table [4] as ORB or, for 
preliminary solutions, ORB:), and in Fig. [5] for the SBs with- 
out orbits. On top of the firm binaries, 5 stars were identified 
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Table 4. Binarity and photometric-variability diagnostics for stars from sample III (as well as binaries or suspected binaries from 
sample IV). 



a. Spectroscopic binaries (marked as black circles in Figs. 2 and 3) 



HD N At Sp o-o (V>) <e> Sb Flag Rem. GCVS a Var. a(Hp) P ph , 







(d) 




(kms -1 ) 


(km s _1 ) 


(kms -1 ) 










(mag) 


83069 


9 


4450 


M2III 


1.63 


0.22 


2.68 


ORB 




NSV 04545 




0.015 


89758 


7 


6632 


M0III 


3.57 


0.25 


3.68 


ORB 


b 


NSV 04829 


SR 


0.019 


108907 


15 


6821 


M4III 


2.57 


0.25 


3.84 


ORB 




CQDra 


cv 


0.053 


111129 


10 


3633 


M2III 


6.81 


0.23 


4.20 


ORB 




BY CVn * 


Irr 


0.070 


113406 


12 


4324 


Ml III 


2.54 


0.23 


2.94 


ORB 




* 




0.017 


132813 


5 


3639 


M4.5III 


4.93 


0.27 


4.89 


ORB 


b 


RRUMi 


SR 


0.107 


137853 


11 


3794 


Ml III 


1.66 


0.22 


2.94 


ORB 




* 




0.014 


159608 


16 


3981 


M2III 


7.87 


0.25 


3.55 


ORB 




* 




0.018 


165374 


37 


7349 


M2III 


5.03 


0.33 


3.61 


ORB 




V980 Her * 


Irr 


0.031 


182190 


24 


4197 


Ml III 


3.15 


0.22 


2.88 


ORB 




* 




0.010 


187372 


16 


5787 


Ml III 


2.42 


0.26 


3.21 


ORB 




* 




0.015 


199871 


10 


3964 


M0III 


7.47 


0.24 


2.33 


ORB 








0.012 


203525 


10 


4741 


M0III 


1.68 


0.23 


0.88 


ORB : 




* 




0.010 


212009 


9 


5008 


M0III 


4.72 


0.23 


3.66 


ORB 




KT Aqr * 




0.028 


212047 


12 


3984 


M4III 


7.22 


0.22 


4.19 


ORB 




PT Aqr 


SR 


0.046 


220088 


15 


4743 


M0III 


1.88 


0.22 


1.67 


ORB 




* 




0.010 


4301 


5 


4607 


M0III 


1.18 


0.25 


2.02 


SB 




* 




0.015 


16058 


5 


4443 


M3III 


3.70 


0.26 


3.65 


SB 




NSV 866 


Lb: 


0.118 


104216 


5 


3626 


M2III 


1.99 


0.27 


3.43 


SB 




FR Cam 


L 


0.064 


134627 


5 


3922 


M0 


1.82 


0.28 


3.31 


SB 




FF Boo * 




0.031 


142804 


5 


4964 


Ml III 


2.64 


0.27 


2.68 


SB 




NSV 7351 


Irr 


0.026 


219215 


7 


4741 


M2III 


0.87 


0.24 


1.72 


SB 


c 


* 




0.010 



a An asterisk in this column means that the variability of the star has been discovered by Hipparcos. 
b HD 89758, HD 132813: an orbit is also available in the literature, as listed in Table[8] 
c HD 219215: Preliminary orbital period: 2500 d 



whose radial-velocity variations are very likely orbital (listed in 
Tabled and displayed in Fig.[6]i. 



3.1.2. Photometric variability 

Although Mira stars have been excluded from our sample, 
some semi-regular variables are nevertheless present. Both 
classes of variables often exhibit pseudo-orbital variations 
cause d by shock wav e s associated with the envelope pulsa- 
tion (lUdrv et all 1 19981 lHatzes & Cochranl Il998t IWoodl 12000; 
Hinkle et al] 120021: ISetiawan et all 120041: (Perekas et all [2006; 
Soszv riski 2007; He kker et al.ll2008l) . For semi -regular variables, 



Hinkl e et al (120021) obtain semi-amplitudes K between 1.6 and 
3.1 km s _1 , whereas for Miras the se mi-amplitudes may r each 
20 km s _1 in the most extreme cases (lAlvarez et al]l200ll) . In 
terms of standard deviations, these values become 1.1, 1.7, and 
14 km s _1 , respectively (remember that for sinusoidal varia- 
tions, cr = Kj \/2). Semi-regular and Mira variables may thus 
be expected to exhibit a^(Vr) values anywhere in the range of 
1 to about 14 km s _1 . 

Moreover, the radial-velocity curves of these stars may often 
be fitt ed by a Keplerian orbit with a period of a few hundred 



be fitt ed by a Keplerian orbit with a period ot a tew hu ndred 
davs (lUdrv" ;t al.ll 1998b iHinkle et alJ l200l IWood et all l2004t 
lLebzelteretaiT l2005). However, it happens that this Keplerian 
solution fits the data for 4 to 10 cycles, and then becomes in- 
valid. Two illuminating examples of this kind of behaviour are 
discussed in the Appendix for stars that do not belong to the 
samples considered in this paper (the Mira S star S UMa and the 



semiregular carbon star X Cnc). If the time span of the radial- 
velocity monitoring is not long enough (i.e., shorter than the 4 
to 10 cycles mentioned above), the inadequacy of the Keplerian 
solution may not be noticed, thus leading to the erroneous suspi- 
cion of binarity. 

It is therefore very important to check that the Keplerian so- 
lutions proposed in Sect. 4 hereafter are valid, by eliminating the 
possibility that they have an intrinsic origin, as often observed 
for Mira and semi-regular variables. For this purpose, we col- 
lected from the literature the photometric properties of all the 
stars listed in Table [4] (including those flagged as having 'radial- 
velocity jitter', among which we should expect to find a large 
fraction of photometric variables). We searched the Hipparcos 
Variability Annex, which has the advantage of being an unbiased 
information source containing all our stars. The results are plot- 
ted in Fig-\3\[cr(Hp)-cro(Vr) diagram], and listed in Table|4] As 
we can see, all stars flagged as spectroscopic binaries (Fig. |3), 
have much larger radial-velocity standard deviations than ex- 
pected from the relation between radial velocity and photomet- 
ric variability for single stars (as reported by Hinkle et al. 1997). 
We also see that this relation exhibits some scatter (Pig. [3j, and 
that stars flagged as 'SB/jitter?' define a rough dividing curve 
between intrinsic and extrinsic radial velocity variations in the 
a(Hp)-ao(Vr) diagram. 

In the column 'Var' of Table |4] we see that almost all stars 
later than M3 are semi-regular (SR) variables, as expected. It 
is very clear that the fraction of SR variables with cr(Hp) > 
0.05 mag increases among stars with (or suspected of having) 
radial- velocity jitter. A specific case (HD 114961, marked as a 
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b. Suspected binaries (marked as black triangles in Figs. 2 and 3) 



HD 


N 


At 


Sp 


ao(Vr) 


Sb 


Flag 


Rem. a 
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Var. 


a(Hp) 


-Pphot 






(d) 




(kms -1 ) 


(km s" 1 ) 










(mag) 


(d) 


65275 


9 


4442 


MO 


0.84 


3.50 


SB? 


III 






0.022 




111946 


9 


5040 


M1/2III 


1.59 


3.02 


SB? 


III 


VW Crv* 


Irr 


0.132 


45.6: 


112142 


5 


4613 


M3III 


1.14 


3.42 


SB? 


IV 


^Vir 


Lb 


0.025 




114961 


5 


4608 


M7III 


3.50 


8.60 


SB? 


IV 


SW Vir 


SRb 


0.198 


150 


189063 


8 


3642 


MOIII 


0.80 


3.16 


SB? 


III d 


* 




0.013 





c. Binarity doubtful (marked as crosses in Figs. 2 and 3) 
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0.77 
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SB/jitter? 


Ill 


AK Cet 


Lb 


0.061 


12292 


5 
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1.00 
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IV 
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17491 


6 


4075 
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1.62 
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IV 
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SRb 


0.157 


24693 


6 


4762 


Ml III 


0.79 


3.11 


ORB/jitter? 


IV b 


GOEri 




0.035 


66875 


8 


4034 


M3III 


1.06 


4.47 


ORB/jitter? 


III b 


BL Cnc 


L 


0.062 


75157 


4 


3626 


M0 


1.01 


3.86 


SB/jitter? 


IV 


FW Cnc* 


Irr 


0.040 


97876 


9 


4841 


M4III 


1.56 


4.22 


ORB/jitter? 


III 


UU Crt* 


SR 


0.115 


129902 


5 


4605 


Ml III 


0.57 


1.87 


SB/jitter? 


IV c 


* 




0.022 


149165 


5 


5087 


M1III 


0.67 


1.82 


SB/jitter? 


IV c 


* 




0.018 


153698 


12 


4209 


M4III 


0.86 


4.09 


SB + jitter? 


III 


* 




0.033 



d. Non-SB (radial-velocity jitter) (marked as open squares in Figs. 2 and 3) 



2411 


6 


4171 


M3III 


1.03 


4.58 


TVPsc 


SR 


0.085 


49.1 


6262 


9 


3702 


M3III 


1.04 


4.51 


V360 And * 


SR 


0.066 




69229 


7 


4819 


Ml III 


0.44 
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0.018 
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8 
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1.04 


5.78 


INHya 


SRb 
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107397 


9 


3746 


M3III 


1.35 


5.42 


RYUMa 
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0.284 


295.7 


112264 


6 


3022 


M5III 


1.08 


6.02 


TUCVn 
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50.0 


115322 


9 


3627 


M6III 


1.49 


5.93 


FH Vir 


SR 
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70.0 


124304 


9 


4361 


M3III 


0.90 


5.05 


EV Vir 


SR 
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120.0 


125357 


8 


4361 


M2/3III 


1.78 


4.35 


MZ Vir 


Irr 


0.151 




166253 


10 


3631 


M0 


1.11 


3.96 


V566 Her 
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0.091 


137.0 


186776 


11 


5981 


M4III 


0.94 


4.47 


V973 Cyg 


SRb 


0.100 


40.0 



a This column specifies whether the binarity suspicion resulted from campaign III or IV (see Table[l} 

6 HD 24693 and HD 66875: Orbital solutions with an eccentricity of about 0.5 and periods of 845 ± 7 and 668 ± 4 d, respectively, are possible, but 
their reality is questionable. The stars are irregular (Lb) variables (varying from Hp — 7. 13 to 7.26 and 5.94 to 6. 1 1 , respectively, in the Hipparcos 
catalogue) falling right on the dividing line in Fig. [2] 

c HD 129902 and HD 149165 are probably spectroscopic binaries, given their location in Fig.[2]just above the dashed line. Unfortunately, the 
radial-velocity data presented in Fig.|6]are too scarce to confirm that conclusion. 
d HD 189063: Orbital period larger than 3600 d 



suspected binary 'SB?' in Tabled} will be discussed in Sect. 13. 31 
Such stars are not very numerous among the ones flagged as 
SB. Therefore, one can be confident that all the stars flagged as 
SB in Tabled are indeed binaries. Only two among those are 
semi-regular variables with large photometric amplitudes (i.e., 
HD 16058 = 15 Tri and HD 132813 = RR UMi), and there are 
good arguments in favour of their binary nature. HD 16058 is a 
newly flagged spectroscopic binary, but it is an X-ray source. 
As we show in Sect. 13.21 X-ray detection in M giants is a 
strong indication of binarity. The binary nature of HP 13 2813 
has already been c laimed by iDettmar & Giesekind d 1983b and 
Batte n & Fletcherl (119861) . who obtained consistent orbital pa- 
rameters (Table[8]l. Our recent radial-velocity measurements are 
consistent with these orbital parameters (see Fig. |9j, thus clearly 
demonstrating the binary nature of HD 132813 from the stability 
of the Keplerian solution. 



3.2. X-ray emission as binarity diagnostic 

X-rays are normally not observed in single M giants lying 
to the right of the so-called dividing line, an almost vertical 
boundary in the HR diagram separating stars with hot coro- 
nae emitting X-rays (to the left) from stars with high mass 
loss (to the right), which are normally not X-ray emitters 
dHiinsch et al.l 1 1998h . Since X-rays may be generated by sev- 
eral physical processes in a binary system, their detection in 
an M giant is a strong argument in favour of its binary na- 
ture. First, X-rays may be produced at the shocks resulting from 
the collision of streams in the complex flow pattern associ- 
ated with wind ac cretion in a detached binary system involv- 
ing an AGB star dTheuns & Jorissenlll993l: iTheuns et al.l l 1996; 
iMastrodemos & Morrislll998l:lJahanara et al.ll2005l) . Second. X- 
rays are generated when the gravitational energy of the M giant 
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wind falling in the potential well of the companion star is con- 
verted into radiative energy when hitting the stellar surface. The 
accretion luminosity will either be radiated away in the form of 
hard X-rays if the infalling matter is optically thin; or if that mat- 
ter is optically thick, half will be converted into thermal energy 
and half will be radiated away in the form of blackbody radiation 
(according to the virial theorem). 

A search for X-ray sources among M giants is therefore 
an efficient way of finding binary systems. The ROSAT all-sky 
survey of X-ray sources detected only 1 1 out of 482 M giants 
of luminosity classe s I to III from the Bright Star Catalogue 
(Huns chet al.| [l998). They are listed in Table [5] along with a 
comment regarding the binary nature. These few detections jus- 
tify using radial-velocity variations to detect binaries among M 
giants in our paper. 

It may be seen from Table [5] that, in most cases, the M gi- 
ants detected as X-ray sources are flagged as binaries (4 from 
the present work, 1 from the literature, and 1 symbiotic). In 
three cases (HD 62898, HD 132813, and HD 155035), the X-ray 
source is offset by more than 1 ' from the optical position of the 
giant, whic h casts doubts on the M giant being the source of the 
X-rays (see [Hiinsch et al.ll 19981 for details). In two other cases, 
there is a visual G-type companion, where X-rays may arise 
from coronal emission. Thus only HD 130144 and HD 150450 
have no satisfactory explanation for the origin of the X-rays. In 
Fig. |3] it is HD 130144 that appears amidst the non-binary M 
giants. 



Fig. 5. Radial-velocity data points for the spectroscopic binaries 
with no orbit available yet. 
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Fig. 6. Radial-velocity data points for the suspected spectro- 
scopic binaries. 



3.3. Special cases 
3.3.1. HD 114961 

Since semi-regular and Mira variables may be expected to ex- 
hibit ero(Vr) values anywhere in the range from 1 to about 
14 km s _1 , the value <Jo(Vr) = 3.1 km s _1 observed for the 
semi -regular variable HD 114961 (SW Vir) in sample II is not 
incompatible with intrinsic variations (although we flagged it as 
a suspected binary; indeed, HD 114961 is the triangle in Fig. [2] 
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Fig. 7. Radial-velocity variations for stars falling close to the di- 
viding line between jitter and orbital motion in Fig. [2] for which 
the nature of the the radial-velocity variations is unclear. The 
vertical scale spans 10 km s _1 in all cases. They are denoted by 
crosses in Fig. [2] 

Table 5. Properties of the 1 1 M giants detec ted as X-ray source s 
by the ROSAT all-sky survey (following iHiinsch et all 1 19981) . 
where for three stars, the X-ray source is offset by more than 
1 ' from the M giant optical position. 
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Fig. 8. Examples of radial- velocity variations for stars (referred 
to by their HD number, as listed in each panel) falling close to 
the dividing line between jitter and orbital motion in Fig. [2] with 
a radial-velocity jitter most likely of intrinsic nature. They are 
denoted by open squares in Figs.|2]and|3] 



with the largest ao(Vr) for Sb = 8.5 km s _1 ). Let us note, 
however, that the radial-velocity data available for this star are 
unfortunately too scarce to infer any periodicity, so no firm con- 
clusion can be reached at this point regarding the intrinsic or or- 
bital origin of the large radial-velocity scatter exhibited by SW 
Vir. 
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Fig. 9. Orbit from iBatten & Fletcherl dl986h for HD 132813. 
Overplotted are our new data points, consistent with this orbit 
from the literature, obtained about two decades ago. 



3.3.2. HD 115521 

HD 115521 is an interesting case, similar to those discussed in 
the appendix. It belongs to sample I (and therefore does not ap- 
pear in Table [4] since it has only been measured twice in our 
observing cam paign), but being a radial-velocity standard star 
till the 1990s (lUdrvet all fl999a). it was measured very fre- 
quently and turned out in the first place to be variable with a 
small amplitude (K ~ 1.5 k m s _1 ) and a period around 500 d 
(Duq uennov & Mavorlll99l1) . The 127 measurements are listed 
in Table[6] only available electronically at the CDS, Strasbourg. 
Later on, variations on a much longer time scale became ap- 
parent, exceeding the measurements' time span of 6358 d. The 
orbital period therefore cannot be determined with good accu- 
racy. Moreover, owing to this insufficient sampling of the or- 
bital cycle, the value adopted for the orbital period strongly in- 
fluences the value derived for the eccentricity. Table [7] lists the 
(pseudo-)orbital elements used to draw the radial-velocity curve 
of Fig. [10] The high uncertainties on the elements of the long- 
period orbit should serve as a reminder that these elements are 
very uncertain. Only lower bounds to the orbital period and the 
eccentricity are therefore listed in Table [8] 

The properties of the short-period variations would imply a 
rather low mass for the companion: assuming a mass of 1.3 M Q 
for the giant, the minimum mass for the companion is 0.054 M , 
corresponding to a brown dwarf. It is not entirely clear whether 
the short-period variations are indeed due to an orbital motion, 
for several reasons. First, as for S UMa and X Cnc discussed in 
the Appendix, the data deviate from the Keplerian solution after 
a dozen cycles, and this is very clearly seen in Fig. [TUJ where 
the first data points do not fall on the solution defined by later 
measurements. Second, with log g = 0.70 and K — 1.2 km s _1 
(Table [7]), the short-period variations of HD 115521 fall on the 
log g — K (K being the semi-amplit ude of the rad i al-vel ocity 
variations) relationship as defined by He kker et alj d2008l) for 




4.4xl0 4 4.6x10" 4.8x10" 5x10" 



JD 

Fig. 10. The radial-velocity measurements of the former radial- 
velocity standard HD 115521. Short-term variations with a pe- 
riod of 475 d are superimposed on variations with a much longer 
time scale. The solid curve is obtained by adopting P = 29208 d 
and e = 0.87 for the long -period orbit. 

K giantfl- The existence of this relation between the ampli- 
tude of the radial-velocity variations and an intrinsic property 
of stars like the surface gravity hints at an intrinsic origin of 
su ch radial- velocity va riations. The relation log g — K found 
by [He kker et al. ( 2008) among K giants (initially suggested by 
lHatzes & CochrarJl 998) continues in the domain of M giants, as 
we discuss in Paper II, where it is shown that the radial-velocity 
standard deviation correlates with the CORAVEL parameter Sb 
(measuring the spectral line width), which is in turn a good mea- 
sure of the stellar radius (see Figs. 1 and 3 of Paper II). 

However, if the shorter-period radial-velocity variations are 
intrinsic, then one expects photometric variations with a pe- 
riod of 475 d; unfortunately, we could not find any mention of 
them, in either the Catalogue of Suspected Variable Stars, where 
HD 1 15521 is entry NSV 6173, or in the Hipparcos Photometry 
Anne x, or in the ASAS Catalogue of Variable Stars dPojmanskil 
l2002h . 

4. Orbital elements of newly-discovered binaries 

The complete set of orbital elements for the 12 newly discov- 
ered binaries are listed in Table [8] Figure [TTJpresents the phase 
diagrams for those firm orbital solutions. The second part of 
Table [8] provides a preliminary period and eccentricity for the 
binary with not enough data points to derive meaningful or- 

5 The gravity estimate is obtained as follows. From V = 4.80 
and the 2MASS value K = 0.47, one derives T eff = 3680 K 
and B C K = 2.74, from the V - K calibrations of iBessell & Wood 
fl984h andlBesselletal.1 fl99g) . The maximum-likelihood distance of 
iFamaev et alj d2005h (168 ± 20 pc) then yields M hol = -2.9, and a 
radius of S4Rq from the Stefan-Boltzmann relationship between lumi- 
nosity, T e ff and radius. Adopting a mass of 1 .3 Mq then yields a gravity 
oflogp = 0.70. 
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Table 6. The 127 radial-velocity measurements for HD 115521. 
The table, only available at the CDS, Strasbourg, lists the Julian 
Date, the radial velocity, and the corresponding error. 



Table 7. Orbital elements for the system HD 115521. It is not 
certain, however, that the short-period variations (A+a) are due 
to an orbital motion (see text). The long period orbit is not well- 
constrained. 





A+a 


Aa+b 




-26.93 ±2.1 





P(d) 


475.1 ± 1.4 


29000 : 


To 


2 449 411 ± 29 


2 419 767 ± 86282 


e 


0.148 ± 0.062 


0.87 ±0.24 


w(°) 


-27.1 ± 22.9 


191.8 ±3.5 




1.25 ±0.07 


4.2 ± 0.5 


f(M) (M@) 


9.3 10" 5 


0.027 


ai sini (AU) 


0.053 


5.57 



(j(O-C) =0.50 km s- 1 



bital solutions. For the sake of completeness, the last part of 
Table [8] collects periods, eccentricities, and mass functions for 
(non-symbiotic) M giants available in the literature, or kindly 
communicated by R. Griffin. This combined data set will be used 
in Paper III to discuss general properties (like the eccentricity- 
period diagram) of systems involving M giant primaries. It must 
be stressed that Table[8]includes orbits neither for symbiotics nor 
for VV-Cephei-like systems (VV Cep, AZ Cas, etc.). A list of or- 
bital e lemen ts for the form e r may b e found in iBelczvnski et all 
d2000l) and [Mikolaiewskal (120031) . iMikoIaiewskal d2007l) and 
iFekel et ail d2007l) provide references for the most recent symbi- 
otic orbits. 

For HD 108907 (4 Dra = CQ Dra), Table lists three 
orbital solutions. The first entry is obtained with our own 
CORAVEL and EL ODIE data alon e . The second entry is a so- 
lution computed by iReimers et al. (1988), based on their 57 
recent data points plus 14 much older ones. The third orbit 
has been computed by merging our 15 data poi nts with the 57 
Camb ridge/CORAVEL/Victoria data points from lReimers et ail 
(119881) . and this combined solution is p resented i n Fig . [121 
This system is of special interest, since iReimersI (1 19851) and 
IReimers et all d 19881) argued that the companion of the red gi- 
ant is a cataclysmic variable, because International Ultraviolet 
Explorer spectra revealed a steep rise shortward of 140 nm along 
with broad lines (full widths of 1000 km s _1 ) of highly excited 
species like He II and C IV. That co nclusion has been cha l lenged, 
howev er, by more recent studies (Wh eatlev et al J 120031: ISkopall 
2005a b) based on ROSAT X-ray observations. They conclude 
instead that the companion is a single WD accreting from the 
wind of its red giant companion, as in normal symbiotic systems. 
The residual radial-velocity jitter of 0.6 km s _1 (listed as cr(0- 
C) in Table[8]l appears normal for a star with Sb = 3.8 km s~\ 
as seen from Fig. [2] 
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Table 8. An extensive list of orbital elements for (non-symbiotic) M giants. 
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Appendix A: Pseudo-orbits among Mira and 
semiregular variables 

This Appendix presents two cases of (suspected) pseudo-orbital variations exhib- 
ited by Mira and semi-regular variables. These stars do not belong to the samples 
studied earlier in this paper. 

A.1. HD 110813 

HD 1 10813 (S UMa) is a Mira S star with a light cycle of 225.9 d (as listed in 
the General Catalogue of Variable Stars - GCVS; the analysis below shows that 
a period of 222 d seems more appro priate), and was considered a spectroscopic 
binary of period 593 ±62 d by[Ulry_et_alj 1 1998). It deserves a follow-up discus- 
sion here, since recent datapoints no longer support this orbital solution; in fact, 



no satisfactory fit can be found to the radial-velocity data. This star provides a 
good illustration of the difficulties encountered while trying to find spectroscopic 
binaries among long-period variable stars. 

Figure I A. 1 1 shows the radial-velocity curve of S UMa and a (pseudo- 
orbital solution based on 17 measurements (from 1987.946 to 2002.464 or JD 
2447 141.743 to JD 2452 444.403; Tables |A~T1 and lA~2l This pseudo-orbit is 
slightly different from the one found by Udrv et al. 1 1998). Although the radial- 
velocity data could be fitted with a period of 576 d for 10 cycles, the last two 
measurements (JD 2452 824.356 and JD 2453 048.562) deviate markedly from 
this solution. The radial-velocity variations cannot therefore be ascribed to or- 
bital motion. In any case, a system with such a short orbital period cannot be 
detached. Using the pe riod-radius rela tionship for Mira stars pulsating in the 
fundamental mode (van Leeuwen et al. 1997), a radius of 258 Rq is inferred 
from the 222-d period, assuming a mass of 1 .5 Mq . Adopting a mass of 1 Mq 
for the companion, for the Roche radius to be larger than the stellar radius re- 
quires an orbital period of at least 1 130 d, which is inconsistent with the observed 
value of 576 d. 

What then is the origin of this 576 d period? Th e Stellingwerf (phase- 
dispersion minimisation) 9 statistics I Stellingwerf 1978) are shown in Fig. IA.2l 
constructed from all datapoints but the last two outliers. It shows that all the 
prominent peaks are combinations of the Mira frequency fx = 1/222 d -1 and 
of the yearly frequency /2 = 1/365.25 d — 1 , or harmonics of f\. In particular 
the 576 d period may be identified with the frequency /i — /2. This finding, along 
with the previous result for the minimum period allowed by the Roche radius, 
definit ely denies the reality of the binary nature of HD 1 10813. Alv arez et alj 
( 200 1 ) have shown that S UMa exhibits an asymmetric cross-correlation dip, as 
usual among Mira variables. Asymmetric profiles observed in Mira variables are 
often associated w ith radial- ve locity variations, which indeed mimick an orbital 
motion (see also Hinkle et al. 2002). 



A.2. HD 76221 

The semiregular c arbon star HD 76221 (X Cnc) behaves similarly to HD 1 10813 
(Udry et al. 1998): a satisfactory (pseudo-)orbital solution with a period of 530 d 
(Table IA.lt could be found with the first 1 3 data points, but is not confirmed 
by the last two data points (bot tom panel of Fig. IA.lt . The phase-dispersion 
minimisation statistics (Fig. |A.3| based on the first 13 data points) reveals several 
harmonics of the 195 d photometric period, but this time, the 530-d periodicity of 
the radial velocities does not seem to be one of these. This 530-d-signal may be 
yet another e xample o f the long secondary periods (of unknown origin) found by 
Houk 1 1963) (see also Hinkle et al. 2002) among SR variables, since the ranges 
of periods, mass functions, and semi-amplitudes found by Hinkle et al. (2002) 
(and listed in Table lA. ll all match the values for HD 7622 1 . Only the eccentricity 
does not conform to the Hinkle et al. range (0.32 - 0.37, with one case at 0.08). 
Soszyhski (2007) also notes that the long secondary periods sometimes undergo 
phase shifts; however, we do not find it very plausible that these photometric and 
spectroscopic phase shifts h ave the same phys ical origin. The phase shifts in the 
light-curve are attributed by Soszyhski 1 2007) to dust clouds in the vicinity of a 



12 

Table 8. (Continued.) 
Preliminary orbit 
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Prieur et al. (2006) 
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5 Sae: A C Aur/VV CeD-like svstem (Reimers & Schroeder 1983) 
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(6) 
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0.032 


Griffin (1990) 
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Remarks: (1) this wor k (2) Rei mers et aD dl988l) (3) combined solution (see text) (4) SS Lep: M4III + accreting Al star JWeltv & Wadell 19951) . 
Verh oelst et alj d2007l) shows, from interferometric measurements of the M-giant radius, t hat th e giant fills its Roche lobe (see also Paper III) (5) 
Possibly a triple system, formerly a radial-velocity standard star; Duquennoy & Mavor ( 1991) found a preliminary 509 d period for the inner 
pair, but the orbital nature of the short-period variations is questionable (Sect. [3~3~2l > (6) HD 190658 = V1472 Aql = HIP 98954: semiregular 
or most likely eclipsin g or ellipsoidal variable with P = 100.37 d, almost exactly half the orbital period ( Samus 1997). Orbital elements from 
Lucke & Mavor ( 1982). A (physical?) companion is present at 2.3" (Tokovinin 1997; Frankowski et al. 2007) 



companion, whereas the phase shifts in radial velocity preclude the orbital nature 
of the radial velocity variations. 
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Table A.l. Pseudo-orbital solutions for the Mira S star HD 1 10 813 (S UMa) and the semi-regular carbon star HD 76221 (X Cnc), 
based on the complete data set but the last two points (see Table |A.2l . 
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a The ranges found by Hinkle et al. ( 2002) for semi-regular variables with long secondary periods. 



Table A.2. Radial velocities for HD 76221 and HD 1 10813. 
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Fig. 12. Top panel: Radial-velocity curve for HD 108907 
(4 Dra), an M giant with a hot compact companion, merging 57 
data points from the Cam bridge, Victoria, and CORAVEL spec- 
trovelocimeters listed in iReimers et alj d 1 9881) (open symbols) 
and the 15 data points from this paper (filled symbols). Bottom 
panel: Phase diagram for the orbital solution merging these two 
data sets. 
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Fig. 11. Orbital phase diagrams for the orbital solutions, accord- 
ing to the solutions listed in Table [8] 
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Fig.A.l. Top panel: Radial-velocity curve for HD 110813, a 
Mira S star with a pulsation period of 225.9 d (according to 
GCVS) and a pseudo-orbital period of 576 d. Bottom panel: 
Same for the carbon star HD 76221, a semi-regular variable with 
a pulsation period of about 195 d (according to GCVS), and 
a pseudo-orbital period of 530 d. For both stars, the last mea- 
surements become discrepant with the (pseudo-)orbital solution 
computed from earlier data points. 
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Fig. A.2. Top panel: Stellingwerf's 1 — 9 statistics (with bins of 
length 0.2 phase cycle offset by 0. 1 phase cycle from the previ- 
ous one; Ni, — 5 and N c ~ 2 adopting Stellingwerf's notations), 
for all data points of HD 1 10813 (S UMa) but the last two (out- 
liers). Middle panel: the 1 — 9 statistics for a sinusoidal signal of 
period 1/fi — 222 d (corresponding to the Mira pulsation pe- 
riod) sampled as the data. Bottom panel: the difference between 
the former two, revealing peaks not accounted for by the aliases 
associated with the 222-d Mira period. In the top panel, several 
linear combinations of f\ (the Mira frequency) and fi (the 1 -yr 
frequency) are indicated. The pseudo-orbital period of 576 d is 
close to (A - /2)- 1 = 555 d. 
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Fig.A.3. Top panel: Same as Fig. IA.2I for the carbon star 
HD 76221, for all data points but the last two (outliers). The solid 
vertical line marks the frequency of the pseudo-orbital variation 
and the dot-dashed line the frequency of light variations. The 
harmonics 1/3, 1/2, 3/4, and 3/2 of the pulsation frequency 1/195 
= 5.1 10~ 3 d _1 are cleary identifiable in the middle and top pan- 
els. Because there is a residual peak remaining in the lower panel 
at the pseudo-orbital frequency of 1/530 = 1.9 10~ 3 d _1 , this 
frequency does not belong to this series of harmonics. 



